HOUSTON  UNI V  TX  DEPT  OF  ELECTRICAL  ENGINEERING  F/G  14/2 

electromagnetic  methods  of  nondestructive  EVALUATION. (U) 

may  61  s  a  LONG.  A  B  EL-KAREH  AFOSR -77-3457 

AFOSR—TR— 82-0674  NL 


UNCLASSIFIED 


AD  All 8587 


m  * 


AEOSR-TR-  8  2-0074 


ELECTROMAGNETIC  METHODS  OF  NONDESTRUCTIVE  EVALUATION 

Final  Report 
through  April  30,  1981 


Sponsored  by  the 

Air  Force  Office  of  Scientific  Research 


a_ 

o 

UJ 


£2 
h  ■- 


under 

Grant  No.  AFOSR  77-3457 


Prepared  by 


Stuart  A.  Long 
A.  B.  El-Kareh 


Department  of  Electrical  Engineering 
University  of  Houston 
Houston,  Texas  77004 


May  18,  1981 


Approved  for  public  release; 
distribution  unlimited. 


Distribution  of  this  document  is  unlimited. 

8°  08  25  002 


'  J  - 


■•Sllfe1-  ?* 


SvMt;  .'.■gT  ,*ff>WMWW^ 


»^ity  classification  of  this  page  fl w>w  PKi^twi^ 

_ REPORT  DOCUMENTATION  PAGE _  I  Before  completSg^orm 

trpnRT  uhm B Um  O  a  m  i  |2.  OOVT  ACCESSION  NO  »•  RECIPIENT'S  CATALOG  NUMBER 

AEOSR-TK-  82-U674  $5*7 


4.  TITLE  (and  Subtltta)  *•  TYPE  OF  REPORT  A  PERIOD  COVEREO 

ELECTROMAGNETIC  METHODS  OF  NONDESTRUCTIVE  EVALUATIC  N  FINAL 

30  Sep  77  -  30  Apr  81 

6.  PERFORMING  O^G.  REPORT  NUMBER 


7.  AUTHORO) 

STUART  A  LONG 
A  B  EL-KAREH 

Y  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  ™" 

UNIVERSITY  OF  HOUSTON 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 

HOUSTON  ,  TX  77004 

It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH/NA 
BOLLING  AFB,  DC  20332 


B.  CONTRACT  OR  GRANT  NUMBERS) 

AFGSR-7 7-3457 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  &  WORK  UNIT  NUMBERS 

61102F 

2307/B2 

12.  REPORT  DATE  ~~~ 

May  18,  1981 

13  NUMBER  OF  PAGES 


14.  MONITORING  AGENCY  NAME  ft  AODRESS^i/  different  from  Controlling  Office;  15.  SECURITY  CLASS,  (of  thia  report) 

UNCLASSIFIED 


\Sm.  DECLASSI  F| CATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thla  Report) 


Approved  for  Public  Release;  Distribution  Unlimited. 


I  17.  DISTRIBUTION  STATEMENT  (of  tha  abatract  entered  In  Block  20,  if  diffarant  from  Report) 


I  18.  SUPPLEMENTARY  NOTES 


It.  KEY  WORDS  (Conlfmi*  on  revaraa  aide  If  nacaaaary  and  Idantlty  by  block  number) 

NONDESTRUCTIVE  EVALUATION 
EDDY  CURRENT  TESTING 
SURFACE  ELECTROMAGNETIC  WAVES 


^  1 20.  ABSTRACT  (Continue  on  rarer  a  a  aide  It  naceaaary  and  Idantlty  by  block  number) 

A  study  of  quantitative  nondestructive  eddy  current  testing  techniques  9tr\ 
r.onferromagnetic  structural  metals  has  been  completed.  In  addition  an  investi¬ 
gation  was  made  of  electromagnetic  surface  wave  propagation  in  a  dielectric 
layer  on  a  conducting  substrate. 


M73 


j6*JN CLASSIFIED _ s _ 

SECURITY  CL ASSl F|C ATlON  OF  THIS  PAjlf1rh#"  Dmt*  **»•*•<> 


TABLE  OF  CONTENTS 


PAGE 


I.  INTRODUCTION .  * 

II.  SUMMARY  OF  PROGRESS .  5 

A.  Eddy  Current  Testing . 5 

B.  Microwave  Testing  .  6 

III.  PUBLICATIONS  AND  PRESENTATIONS  .  7 

IV.  SCIENTIFIC  PERSONNEL  .  8 

Appendix  A  -  Impedance  of  a  Loop  Surrounding  a  Conducting  Cylinder  .  9 

Appendix  B  -  The  Impedance  of  a  Single-Turn  Coil  Near  a  Conducting 

Half  Space .  15 

Appendix  C  -  Nondestructive  Evaluation  of  Dielectric  Layers  on  Con¬ 
ductive  Substrates  by  Microwave  Surface  Electromagnetic 
Waves .  23 

Appendix  D  -  The  Change  In  Impedance  of  a  Single-Turn  Coil  Due  to  a 

Void  in  a  Conducting  Half  Space .  30 


AIR  FORCE  O'FFICE  of  scientific  RESEARCH  { AFSC  ' 

NOTICE  OF  TRANSMITTAL  TO  DTIC 

This  technical  h-s  *  r*vi**«d  »nd  is 


approved  ? 

Distribut Ic 


•  '  ”  ’  •  *  ^  IAA  AT  R  19  3-12. 

uni imited. 


MATTHEW  J.  KERPER 

Chief,  Technical  Inf ormation  Division 


I. 


Introduction 


The  general  objective  of  this  program  was  to  further  the  development 
of  nondestructive  evaluation,  particularly  in  relation  to  its  use  as  a 
means  of  quantitatively  characterizing  performance  related  properties  of 
structural  materials.  The  principal  two  areas  of  emphasis  have  been  low 
frequency  eddy  current  testing  methods  for  nonferromagnetic  metals  and 
microwave  testing  of  dielectric  layers  on  conducting  substrates  using 
surface  electromagnetic  waves. 
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II.  Summary  of  Progress 

A.  Eddy  Current  Testing 

The  excitation  of  eddy  currents  in  materials  to  detect  flaws  is 
well  developed  in  practice.  The  theoretical  solutions,  however,  of  even 
the  most  basic  geometries,  which  even  remotely  resemble  practical  testing 
situations,  have  not  been  attempted  until  recently.  The  numerical  solutions 
of  Dodd  and  Deeds  [Journal  of  Applied  Physics,  Vol.  30,  pp.  2823-2838,  1968] 
and  the  analytical  work  by  Zaman,  Gardner  and  Long  for  both  cylindrical 
[IEEE  Transactions  on  Instrumentation  and  Measurements,  March  1981]  and 
planar  [Journal  of  Nondestructive  Evaluation,  1981 1  geometries  are  the 
first  real  attempts  to  attack  the  basic  eddy  current  problem  on  a  theo¬ 
retical  level.  The  results  of  these  studies  have  application  in  many 
practical  cases  where  eddy  current  methods  have  been  used  for  years. 

The  case  of  a  single-turn  loop  surrounding  an  imperfectly  con¬ 
ducting  cylinder  has  been  solved  for  a  slightly  restrictive  set  of  physical 
parameters.  The  change  In  complex  impedance  of  the  coil  was  calculated  as 
a  function  of  the  geometry  of  the  problem  (radii  of  the  coil  and  core)  and 
of  the  material  properties  of  the  core  (conductivity) .  (See  Appendix  A) 

In  a  similar  fashion  the  impedance  of  a  loop  parallel  to  and  near 
an  infinitely  large  planar  conductor  was  calculated.  This  change  in  com¬ 
plex  impedance  was  found  as  a  function  of  the  size  of  the  coil,  the  lift¬ 
off  distance  and  the  conductivity  of  the  material.  Again  these  results 
bear  direct  application  for  practical  testing  situations  employing  planar 
geometries.  (See  Appendix  B) 

The  results  of  these  previous  investigations  may  also  be  used  to 
calculate  the  change  in  impedance  due  to  a  flaw  in  the  conducting  material. 


A  first  approximation  using  only  the  fields  in  the  unflawed  sample  has  been 
developed  for  the  usual  case  of  a  single  coil  eddy  current  system  in  the 
presence  of  a  conducting  half  space.  (See  Appendix  D  ) 

B.  Microwave  Testing 

The  investigation  of  the  use  of  surface  electromagnetic  waves  to 
measure  the  thickness  and  dielectric  constant  of  a  dielectric  layer  sup¬ 
ported  by  a  planar  conductive  substrate  has  been  completed.  The  theoretical 
results  show  that  the  thickness  and  the  dielectric  constant  can  be  measured 
independently  by  exciting  a  surface  electromagnetic  wave  along  the  layer 
and  the  subsequent  measurements  of  its  propagation  characteristics*  The 
results  of  this  theory  have  been  tested  by  an  experimental  investigation 
of  surface  waves  traveling  along  1-2  cm  thicknesses  of  layers  of  dielectric 
materials  at  a  frequency  of  10  GHz.  The  resulting  predictions  of  thickness 
and  dielectric  constant  are  found  to  be  quite  accurate  when  applied  to 
samples  of  known  physical  properties.  The  measurement  of  thinner  layers 
may  be  accomplished  by  simply  increasing  the  frequency  of  operation. 

Some  addition  complexity  of  the  actual  experimental  apparatus 
will  result,  but  the  same  theory  will  still  be  valid.  Technical  details 
of  both  the  theoretical  work  and  the  experimental  set-up  including  a  com¬ 
parison  of  theory  and  experiment  is  included  as  Appendix  C,  A  more  detailed 
paper  is  presently  under  preparation  and  will  be  submitted  to  a  suitable 
journal  for  publication. 
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Impedance  of  a  Loop  Surrounding  a  Conducting 

Cylinder 

AFROZ  J  M.  ZAMAN,  STUART  A.  LONG,  member,  IEEE,  and  C.  GERALD  GARDNER 


Abstrnct— The  change  in  complex  impedance  between  an  ideal 
one-turn  coil  surrounding  and  coaxial  with  an  infinitely  long  circular 
cylinder  of  conductivity  a  and  permeability  p  and  a  simitar  coil  without 
the  core  has  been  calculated.  From  the  exact  expression,  a  power  series 
in  the  quantity  6/b  (6  *  skin  depth;  b  ~  radius  of  core)  has  been  de¬ 
veloped.  From  this  result,  the  change  in  impedance  of  a  physically 
realistic  multi  turn  coll  can  be  estimated.  The  theory  permits  a  rational 
approach  to  optimization  of  the  design  of  eddy-current  test  coils  and 
provides  a  possible  means  of  detecting  changes  in  the  radius  and  con¬ 
ductivity  of  the  cylinder. 

I.  Introduction 

A  NOTABLE  OMISSION  to  the  present  body  of  k'nowl- 
.  edge  dealing  with  eddy-current  testing  is  an  adequate 
theoretical  basis  for  the  interpretation  of  changes  in  the  im¬ 
pedance  of  the  test  coil.  This  deficiency  remains,  even  though 
the  fundamental  theory  is  well  established,  owing  to  the 
mathematical  difficulties  involved  in  solving  the  equations  for 
practical  test-coil  and  specimen  configurations.  A  complete 
solution  in  analytical  form  seems  to  exist  for  only  a  few  ide¬ 
alized  cases  which  do  not  necessarily  approximate  practical 
problems  of  current  interest. 

The  problem  selected  for  study  in  this  investigation  is  that 
of  an  idealized  one-turn  coil  (or  loop)  around  and  coaxial  with 
a  long,  solid,  electrically  conducting  cylinder.  This  arrange¬ 
ment  is  illustrated  in  Fig.  1  and  shows  the  loop  with  radius  a 
and  the  core  with  radius  b  and  conductivity  a.  The  theoretical 
treatment  will  assume  that  the  core  is  infinitely  long.  (This 
approximation  should  produce  very  small  errors  as  long  as  the 
distance  from  the  position  of  the  loop  to  either  end  of  the  core 
is  large  compared  to  the  dimensions  of  the  loop  itself.)  This 
problem  has  the  advantage  of  being  simple  enough  to  permit 
a  meaningful  approximate  solution  to  be  found  while  still 
corresponding  to  a  practical  eddy-current  testing  situation.  The 
results  show  how  the  complex  impedance  of  the  test  coil 
changes  when  a  cylindrical  specimen  is  placed  inside  the  loop 
and  how  this  impedance  is  a  function  of  the  geometrical  and 
material  parameters  of  the  cylindrical  core. 

II.  Theory 

A  theoretical  treatment  of  a  geometrically  similar  problem 
has  been  previously  reported  by  Islam  [1J.  In  that  work, 
however,  the  emphasis  was  on  the  radiation  properties  of  the 

Manu*cript  received  October  1, 1979;  revised  September  9, 1980.  This  work 
wav  supported  in  part  by  the  US.  .  Force  Office  of  Scientific  Research 
under  Grant  77-3457 

The  authors  are  with  the  Electrical  Engineering  Department,  University 
of  Houston,  Houston,  TX  77004. 


Fig.  I.  Loop  with  a  cylindrical  core. 


configuration  and  thus  only  a  high-frequency  approximation 
was  attempted  for  the  case  of  a  magnetically  permeable  cy¬ 
lindrical  core.  The  case  of  interest  in  this  work,  that  of  an 
electrically  conducting  core  at  much  lower  frequencies,  may 
be  attacked  in  a  similar  fashion  but  is  essentially  a  completely 
different  problem.  A  previous  paper  by  Dodd  and  Deeds  [2| 
investigates  this  configuration  but  its  solutions  are  left  in  the 
form  of  complicated  integrals  which  must  be  evaluated  nu¬ 
merically.  From  Maxwell's  equations  for  time-harmonic  fields 
one  may  derive  the  wave  equation  for  the  magnetic  vector 
potential  A  in  terms  of  the  impressed  current  density  J. 

V2A  +  *2A  =  J.  (1) 

Using  the  standard  eddy-current  approximation  of  neglecting 
the  displacement  current  terms  and  recognizing  that  the  vector 
potential  has  only  a  <f>  component  which  depends  on  r  and  r, 
the  left-hand  side  of  the  equation  becomes 


Equating  this  expression  to  the  source  terms  due  to  the  loop 
current  and  the  induced  eddy  currents  one  obtains  the  equation 
for  the  vector  potential  in  each  of  the  regions  shown  in  Fig. 
I. 
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U-l*  I  <M*  _  T,  \, 

dr  r  (ir  r  ■  <)z  2 


C  4  =  Mob)" 


- /ut)io<Hz)h(a  -  r), 

fa.' fj  CT.-f  ,. 


where  Mn  is  the  permeability  of  tree  space,  /<>  is  the  magnitude 
of  t he  impressed  loop  current,  u?  is  the  angular  frequency  of 
the  lime-harmonic  fields,  and  ^  and  n  are  the  permeability  and 
conductivity  of  the  core  The  presence  of  the  impressed  loop 
current  at  r  =  0,  r  =  a  is  represented  by  the  two  b  func¬ 
tions 

T he  solution  to  the  equation  max  be  found  using  a  cosine 
transform. 


.[  kbl0(kb)lt(y  )-•)/„(->  )l](kh)  j  . 

4 - - -  A  ) 

kbl\ (>  )k{)(kb)  4  4/n(y  ) A  i ( /ccr ) 


w  here  /oand  A'oarc  the  modified  Bessel  functions  of  order  zero 
and  =  \  (Ah)2  4  /(sh):  Using  this  expression  the  value  of 
along  the  loop  at  r  =  0,  r  =  a  can  be  found 


A+ia.  0)  = 


I  \  (ka)K \{ka)  dk  4 


■s:\ 


[  khi{){  A/»)/|( *> )  —  y  / o( )/ 1 ( kb)  1  <, 

-  - - —  ^  r(  Au)  u  A 

[Ah/ it") ) A n(  A /> )  4  y/0(4  )A  i(Ah)j 


r.  r) 


^  r,  A  )  cos  kz  dk 


with  the  following  functions  defined  for  each  of  l he  three  re- 


r,/iK  (A  •  4  j k  *V|.  r  <  b 

tfl  rk)  =  <  (4/,  Ur)  4  (\K  ,Ur).  b  <  r  <  a  Is) 

C4A  iUrf  r  >  a 

where  / 1  and  k  \  arc  the  modified  Bessel  functions  of  order  one, 
k  2  =  u.'u ^ ,  and  ( *t.  (4.  U  *.  and  (  4  arc  constants  to  be  deter¬ 
mined  In  the  boundarx  conditions 

Since  the  quantii \  of  primary  interest  is  the  vector  potential 
in  the  vicinity  ot  the  loop,  the  simplest  expression  is  that  for 
Region  111  for  whk  nix  (  4  needs  to  be  found  from  the 
standard  boundarx  conditions 


The  first  of  these  integrals  can  be  shown  to  be  exactly  the 
contribution  to  the  vector  potential  due  to  the  loop  itself  if  the 
core  were  not  present  at  a  If  (This  term  is  singular  in  nature.) 
T  he  second  integral  is  the  contribution  due  to  the  eddy  currents 
and  thus  represents  the  difference  in  the  vector  potential  wiib 
and  without  the  core  present  This  term  now  called  A  4^  may 
be  expanded  as  an  asymptotic  series 
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1 1  should  be  noted  that  the  dependence  of  A.T^  on  the  ma- 
tcrial  parameters  /u  and  a  of  the  core  is  completely  contained 
:n  the  /  t  y  )  term  f  or  certain  cases  of  practical  interest,  k  is 
large  compared  to  values  of  A  within  the  range  for  which  the 
i  nice  rands  in  ( S  >  are  appreciable,  lor  such  cases,  7  (and  hence 
/ 1  y  it,  is  nearly  independent  of  A  Indeed,  as  k  becomes  suf- 
f.viently  large  t corresponding  to  a  core  conductivity  ap- 
pioaching  infinity)  we  have 

y  *  x  jkH  =  ( l  +  j){h/b) 

a  here  •>  ~  \  2  u.>ur  ,  the  skin  depth.  We,  therefore,  expect 
that  lor  small  but  nonvanishing  skin  depth,  7(7)  may  be  ad- 
equatek  represented  b\  an  asymptotic  series  in  (b/b).  Using 
the  standard  as-  mptotic  series  for  the  modified  Bessel  func¬ 
tion^  we  tind 


*> 
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1  \p,i  tiding  y  1  as  a  power  series  in  A /k.  and  neglecting  terms 
.  *1  order  lug  her  than  ( a  ■  /> )  \  we  find 
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The  apparent  change  in  the  driving  point  impedance  <>!  ihc  loop 
is 


A/  *  AR  +  fu>Al  = 


jw  A.  Aelita 

U) 


(17) 


From  the  previous  expressions,  the  changes  in  inductance  and 
resistance  can  be  found  to  third  order  in  n/h 


^(f  +  ,,0> 


Note  th.it  the  lowest  order  term  affected  by  the  k  dependence 
-1  y  is  the  1 0  h)x  term.  For  the  range  of  geometrical  and  ma- 
icri.il  parameters  considered  in  subsequent  numerical  work, 
the  cuhk  term  contributes  negligibly  to  the  final  result. 

One  ma\  then  separate  the  contributions  to  A.-f^  into  real 
md  imaginary  parts 


A  l 


Uni  it  it 

»  U 


-An 


O' 


2  \h. 


(IX) 


A/e  =  2w/i°ft  (2  d 


t  pi-’ 

2  A 


1  .  M  ,  .  / « 

2  UJ  ‘ l h 


A\ 


1 19) 


To  obtain  numerical  values  for  A/,  and  A R  it  is  first  necessary 
to  evaluate  the  integrals  A0,  A’|,  /Vs,  S  \%  and  V4  Although 
they  cannot  be  evaluated  analytically  they  can  be  found  nu¬ 
merically  for  fixed  values  of  the  geometrical  parameter  u/b 
Once  these  integrals  are  evaluated,  the  expressions  are  each 
seen  to  be  a  power  scries  in  the  parameter  b/h  which  contains 
the  electrical  properties  of  the  core  material  One  should  note 
that  for  the  case  of  a  perfectly  conducting  core  (i.e  .  hjb  = 
0) 


A/  =  ,\  P 

b  h 


where  th,  N'l  low  mg  integrals  have  been  defined  and  are  seen 
to  be  onk  a  luik  lion  ol  the  ratio  ujb  ( rf  =  kb) 


A  d  >?) 


dr] 


02) 


anti  AR  —  0  This  is  a  reasonable  result  which  shows  a  dci  rease 
in  the  inductance  but  no  change  in  the  resistance  since  no  losses 
are  possible.  T  he  effect  of  a  large  but  finite  conductivity  is  seen 
to  diminish  the  amount  of  decrease  found  for  the  perfectly 
conducting  case  and  to  add  a  finite,  positive  apparent  resis¬ 
tance 


.  :\r\  »rT  n  a 
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III  Rksi  I  TS 

1  v)  facilitate  the  evaluation  of  A/,  and  A/?  for  practical  cases 
the  integrals  \n.  V,.  /Vi.  Af*.  and  /V4  were  evaluated  for  several 
values  of  ajb  varying  from  a  value  of  1 .05  to  2.0.  Using  these 
results,  the  values  of  AA!  and  A L  can  be  calculated  through 
terms  of  order  { A/A)\  Accuracy  of  the  results  depends  criti¬ 
cal  ! i \  on  the  assumption  that  A/A  remains  small  with  respect 
to  one 

f  or  1  lie  case  of  an  aluminum  core  with  a  4- in  diameter  we 
tind  n  *  0.082<>  \  /w  hie h  for/  =  50  kHz  yields  a  skin  depth 
<>  =  0.37  mm  resulting  in  a  value  of  A/A  =  0.0388.  Thus  for  this 
practical  case  we  are  well  within  the  assumptions  used  in  the 
derivations. 

To  generalize  the  results  somewhat,  the  normalized  quan¬ 
tities  AAV  up i/  and  A.V/u have  been  plotted  in  the  re- 
m .’.niing  figures  (A.Y  —  u'A/  ).  It  should  be  noted  that  each  of 
these  quantities  is  unitless.  The  most  obvious  graphical  pre¬ 
sentations  would  be  those  of  A R  and  A/,  versus  the  geomet- 
neal  parameter  ajb  and  the  material  parameter  A/A.  Unfor¬ 
tunately.  this  direct  approach  docs  not  correspond  to  the  usual 
parameters  which  ma>  be  subject  to  change.  Assuming  that 
the  practical  testing  situation  consists  of  a  cylindrical  sample 
mov  mg  through  a  fixed  coil,  the  quantities  which  may  change 
are  actual!)  the  radius  of  the  core  A  and  the  conductivity  of  the 
core  material  rr  lo  illustrate  the  changes  in  impedance  for 
\  a  nations  in  A  about  a  nominal  radius  A(),  the  graphs  in  Figs. 
2  and  *  are  shown  In  Fig  2.  the  normalized  change  in  resis¬ 
tance  is  shown  versus  the  quantity  A/A,»  The  nominal  radius 
A()  may  be  an\  value  as  long  as  our  restriction  of  A/A0  «  I  is 
satisfied  The  vV  *  of  a/bn  ~  1.25  was  chosen  to  be  repre¬ 
sentative  of  a  real  mil  design  which  couples  strongly  with  the 
core  A  family  of  curves  is  also  shown  for  several  values  of  A/A0 


2 

k 


f  ig  5.  Normalized  ren.siancc  versus  conductivity 


°o 


fig  6  Normalized  reactance  versus  conductivity 


It  is  noted  that  all  the  curves  approach  zero  as  A/Ao  is  decreased 
and  become  very  large  as  A/A0  approaches  1.25  which  is  the 
position  of  the  driving  loop.  A  similar  set  of  curves  is  shown  in 
Fig.  3  for  the  change  in  reactance.  Again  as  expected,  the 
change  in  reactance  approaches  zero  as  the  core  radius  de¬ 
creases  and  becomes  a  very  large  negative  value  for  A/An  near 
l  25.  It  should  be  noted  that  changes  in  the  parameter  A/A0 
have  a  relatively  small  effect  on  AA  as  compared  to  their  effect 
on  A R  The  same  functional  dependence  is  also  illustrated  in 
Fig.  4  The  solid  curves  s^ow  the  normalized  resistance  plotted 
versus  the  normalized  reactance  as  A/A0  is  varied.  Changes  in 
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f  ij»  7  Resistance  versus  reactance  as  conductivity  changes. 


the  complex  impedance  can  be  seen  for  varying  radii  for  each 
of  four  values  of  o/Hq. 

The  effect  of  changes  in  conductivity  of  the  core  on  the  re¬ 
sist  a  nee  and  the  reactance  are  shown  in  Figs.  5  and  6.  The 


conductivity  is  again  normalized  with  respect  to  a  <To  near  that 
of  aluminum.  (However,  cr0  is  actually  arbitrary  as  long  as  the 
condition  that  b/b  «  I  is  still  satisfied.)  The  limiting  behavior 
is  again  logical  showing  the  resistance  approaching  zero  for 
large  conductivities  and  zero  for  very  small  values.  The  reac¬ 
tance  is  seen  to  approach  the  “perfect  conductor”  case  as  a 
increases.  The  region  where  a  becomes  small  violates  the  as¬ 
sumption  on  b/b  and,  therefore,  the  behavior  of  these  curves 
has  no  meaning  in  this  region.  The  resistance  versus  the  reac¬ 
tance  is  shown  in  the  solid  lines  of  Fig.  7.  As  the  conductivity 
decreases  from  the  perfect  conductor  case,  AR  is  seen  to  in¬ 
crease  while  AA'  becomes  less  negative.  Each  of  these  curves 
terminates  in  the  region  where  the  assumption  that  b/b  «  I 
begins  to  break  down. 

With  the  aid  of  Figs.  4  and  7  one  may  ascertain  the  behavior 
of  changes  in  both  the  resistance  and  reactance  for  any  percent 
change  in  either  the  radius  of  the  core  of  its  conductivity.  The 
functional  change  in  impedance  is  quite  different  for  the  two 
parameters.  This  characteristic  may,  therefore,  be  utilized  in 
practical  testing  to  determine  changes  in  sample  radius  and 
conductivity  of  rods  and  tubes  w  ith  wall  thickness  greater  than 
several  skin  depths. 
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The  change  in  complex  impedance  between  an  ideal  one-turn  circular  coil  located  above  and 

{parallel  to  a  conducting  half  space  with  respect  to  a  similar  isolated  coil  has  been  calculated 

From  this  result,  a  series  expansion  of  the  integrand  allows  the  solution  to  be  approximated 
by  terms  expressed  as  complete  elliptic  integrals.  Results  have  been  calculated  for  the  change 
in  impedance  as  a  function  of  the  liftoff  distance  and  the  conductivity  of  the  half  space  for  a 
I  coil  of  representative  radius. 


KEY  WORDS:  eddy  current .  impedance;  liftoff;  conductivity 


1.  INTRODUCTION 

The  eddy  current  method  of  nondestructive 
evaluation  entails  the  induction  of  eddy  currents  in  a 
-  conductive  test  object  by  a  time-varying  field  pro 
]  duced  by  a  suitable  distribution  of  impressed  currents 
(via  an  excitation  or  primary  coil),  and  the  detection 
of  the  resultant  field,  usually  by  an  inductive  search 
i  coil,  which  may  be  either  a  separate  secondary  coil  or 

1  the  primary  coil  itself  (see  Fig.  1).  The  method  is 

ordinarily  used  at  frequencies  sufficiently  low  to  ne- 

I  gleet  effects  due  to  displacement  current;  hence  a 
theoretical  analysis  entails  calculating  either  a  trans¬ 
fer  impedance  for  a  primary  coil  and  secondary  coil 

I  in  the  presence  of  the  test  object,  or  the  calculation  of 
the  self-impedance  of  a  primary  coil  in  the  presence 
of  the  test  object.  In  practice,  one  often  needs  only 
the  change  in  impedance  produced  by  the  test  object 

I  or  by  changes  in  the  nominal  properties  of  the  test 
object  (e.g.,  changes  in  its  geometry  or  position  with 
respect  to  the  test  coil  or  coils,  or  distributed  or 

I  localized  changes  in  the  resistivity  of  the  test  object). 
The  most  general  case,  allowing  arbitrary  configura- 

I1  Department  of  Elcctncal  Engineering,  University  of  Houston, 
Houston,  Texas  77004 


lions  of  primary  and  secondary  coils  and  arbitrary 
test  objects,  can  be  handled  only  by  numerical  meth¬ 
ods.  Certain  idealized  arrangements  can  be  treated 
analytically  either  exactly  or  in  useful  approximation. 
In  virtually  all  cases  of  practical  interest,  the  analysis 
eventually  reduces  to  the  evaluation  of  certain  in¬ 
tegrals  that  cannot  be  expressed  in  closed  form  in 
terms  of  standard  transcendental  functions. 

In  this  paper,  we  discuss  the  case  of  a  one-turn 
circular  coil  located  above  and  parallel  to  the  surface 
of  a  homogeneous  conductive  half  space.  From  the 
standard  boundary  value  problem  approach,  we  ob¬ 
tain  the  general  expression  for  the  change  in  coil 
impedance,  A Z,  produced  by  the  half  space;  AZ  is 
given  in  terms  of  an  integral  over  a  separation  param¬ 
eter.  A  series  expansion  of  one  term  in  the  integrand 
permits  the  integral  to  be  expressed  as  a  series  of 
terms,  each  of  which  is  expressible  in  terms  of  com¬ 
plete  elliptic  integrals.  The  leading  terms  of  this  series 
approximate  AZ  asymptotically  for  sufficiently  small 
values  of  skin  depth  of  the  half  space. 

The  problem  addressed  here  has  previously  been 
treated  by  Cheng,0 1  who  evaluated  AZ  by  numerical 
methods  for  various  choices  of  the  relevant  parame¬ 
ters.  Similarly,  Dodd  and  Deeds(2)  have  devised  a 
digital  computer  program  capable  of  handling  cir- 
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cular  test  coils  in  the  presence  of  layered  planar  and 
coaxial  cylindrical  test  objects.  Such  brute  force 
numerical  procedures  are  valuable  for  design  pur¬ 
poses.  but  have  the  disadvantage  of  somewhat  con¬ 
cealing  the  essentially  simple  manner  in  which  the 
final  result  depends  upon  the  parameters  of  the  prob¬ 
lem.  The  approach  taken  here,  while  less  universal 
than  the  purely  numerical  approach,  results  in  rela¬ 
tively  simple,  though  approximate  and  restricted,  for¬ 
mulas  for  1Z  in  terms  of  (he  basic  parameters  of  the 
problem  For  illustrative  and  comparative  purposes, 
some  selected  numerical  examples  are  also  given. 


2.  THEORETICAL  ANALYSIS 


k2 A  term  is  negligible  for  z>0,  we  have 


i 

3  r 


3  z2 


for  z>0 


(2) 


and,  with  k2  —  —  yto/io  for  z<0, 


i  dAt 

- -  -I-  —  - -  -1 - " 

dr2  r  dr  3*2 


A* 

—  — / uiloA  4  =  0  for  z  <  0 


(3) 


Solving  by  the  separation  of  variables  technique 
yields  the  following  expression  for  the  general  solu¬ 
tion  to  Eqs.  (2)  and  (3): 

A„(r.  z)  =  f*[A(a)ea!  +  B(a)e  ~“] 

X  [C(a)7,(ar)-f  Z)(a)Yj(ar)]  da  (4) 


where  a  is  the  separation  constant.  Since  z  may 
become  infinitely  large  in  the  region  z>/,  the  coeffi¬ 
cient  A(a)  must  equal  zero.  Similarly,  in  the  region 
z < 0,  B(a)  must  also  equal  zero,  and  since  the  origin 
is  included  in  all  regions,  then  D(a)  must  equal  zero 
in  each: 

/L,(r,z)=  f  Bxe~azJy{otr)  da 

A> 

z>l>  0  (5) 

r,z)=r[C2e°'+  B2e  "«]  7,(  ar )  da 

A) 

l>z>  0  (6) 


The  basic  geometry  of  the  problem  is  shown  in 
Fig  1.  and  consists  of  a  loop  of  radius  r0  oriented 
parallel  to  and  at  a  distance  /  above  a  homogeneous 
half  space  of  conductivity  a.  Beginning  with  the  basic 
equation  for  the  vector  potential, 

V2A  +  *2A  =  — -S(l-z)  (1) 

and  noting  the  symmetry  of  the  problem,  it  is  seen 
that  the  only  component  of  the  vector  potential  pre¬ 
sent  is  the  circumferential  component,  and  that 
A4  is  a  function  of  r  and  z  only.  Making  the  usual 
low-frequency,  quasi-static  approximation  that  the 


'4*3(r’z):=  f  C^ea'zJy{ar)da 

A> 

z<0  (7) 

where  a2  =a2  +ywjia. 

Since  the  electric  field  is  proportional  to  A+y  the 
boundary  conditions  for  the  tangential  electric  field 
can  be  satisfied  by  equating  the  values  of  A ^  at  the 
z  =  l  plane: 

rB,e-'M ar)Ja  =  /'°0(C!<’o,  +  B2e  ar ) da 

A)  A> 

(8) 
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Multiplying  both  sides  by  the  integral  operator 
/o“{  '  ’  ’  )Jx(a'r)rdr  and  using  the  Fourier-Bessel  iden- 


tity,(3) 


/■OO  ,  Ola"' 

f  Ji(ar)Jl(a'r)rdr= - 

4)  a 

gives  the  following  result: 


a —  a') 


-a7 

a  a 


Bxe~al  C2  ..  B2 
- — —  =  -Teal+-fe~al 


or 


V  a‘  =  C2ea(  +  B2e 


(9) 

(10) 

(11) 


The  two  terms  in  the  square  brackets  represent,  re¬ 
spectively,  the  vector  potential  due  to  the  loop  itself 
and  that  due  to  the  currents  induced  in  the  conduct¬ 
ing  plane.  This  second  term,  due  to  the  conductive 
half  space,  will  produce  the  change  in  impedance 
from  the  case  of  the  isolated  loop  to  the  case  of  the 
loop  near  the  plane.  This  change  in  vector  potential  is 
thus  given  by  this  second  term: 


The  radial  component  of  the  magnetic  field  can  also 
be  found  from  the  vector  potential,  Hr  =  —(8 /<Sz)A^ 
Hr  is  discontinuous  at  the  position  of  the  loop  (r= 
rn,  r  —  /)  by  an  amount  equal  to  the  surface  current 
density  there: 


8  8  ' 
d;  A+'  +  dz  A*2 


■ p/«(r-r0 )  (12) 


or 


This  change  in  vector  potential  can  be  used  to 
calculate  the  change  in  impedance  due  to  the  pres¬ 
ence  of  the  conductor  by  integrating  the  tangential 
electric  field  around  the  position  of  the  loop: 

AZ= - - - -ju— j-&A+x{r0,l)  (19) 

since  —  JvtiA+xfc  then 


Bxe*1  =  C2e" -tiIrzJx(ar0)  (13) 

Similarly,  the  boundary  conditions  may  also  be  ap¬ 
plied  at  ^  =  0.  where  both  E*  and  Hr  are  continuous, 
yielding 


c2+b2=c3 

04) 

and 

C2-S2  =  jC3 

(15) 

These  four  equations  (11,  13,  14,  and  15)  can  then  be 
solved  for  the  constants  /?,,  C2,  B2>  and  C3  and  used 
in  Eqs.  (5),  (6),  and  (7)  to  evaluate  the  vector  poten¬ 
tial. 

Since  our  principle  interest  lies  in  evaluating  the 
vector  potential  at  the  location  of  the  loop,  the  most 
direct  route  is  to  evaluate  the  constant  fl,: 


filr0Jx{ar0) 


ea,+e 


-./0~(ai/tt)) 

(l  +  (a,/a)) 


(16) 


AZ  =  ir«r0Vjf  J*( ar0 )e  2o/(^~)^«  (20) 

The  integrand  factor  (a  — a,)/(a  +  a,),  essentially  a 
reflection  factor,  has  modulus  equal  to  or  less  than 
unity,  the  extreme  value  being  assumed  for  a  =  0  and 
a  —  oo.  The  integrand  factor  J\(ar0)  guarantees  that 
the  value  of  the  integral  is  negligibly  affected  by 
values  of  a  greater  than  about  10/r0.  Practical  values 
of  r0  are  usually  of  the  order  of  10'2m.  For  such 
values  of  r0,  the  importance  range  for  a  is  103 

m\  while  the  quantity  w/ia  [  =  2/(skin  depth)2]  is, 
in  many  practical  cases,  of  the  order  of  107  (e.g.,  for 
aluminum  at  50  KHz,  ioya=  1.5X  107).  For  such 
cases,  a2/w/ia< 0.1,  and  (a~a}  )/(a  +  a,)  may  be 
expanded  as  a  power  series  in  a 

£^  =  _,  +  A£_2ii\  ... 

«  +  fj  «  J  K1 


Thus 


=  -1+0 -j)(a6)+j(ai?+  ■■■  (21) 

where  8  =  y/2/io/ic  ,  and  k  =  Jujio  . 

We  expect  the  series  above  to  converge  rapidly 
provided  a6<fc  1.  As  we  shall  presently  show,  it  is 
convenient  to  adopt  r0  as  a  characteristic  length. 
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Since  the  value  of  A Z  is  determined  almost  entirely 
by  values  of  a  for  which  ar0  <10,  we  have  rapid 
convergence  of  the  integrated  series  if  6/>0<  1/10. 
Separating  A Z  into  real  and  imaginary  parts,  we  have 


\Z  =  bR+jbX 

(22) 

AX=  —  jnoTrr02 j  j  J2(ar0)e~2al  da 

-  j  8J2(  ar0)e~ 2ata  da  J 

(23) 

A/?  =  jicj7rr02jy*  8J2(ar0)e  ~2ata da 

-  J  82J2(ar0)e~2aia2  </aJ 

(24) 

These  changes  in  resistance  and  reactance 
represented  by  three  integrals: 

can  be 

AA=  —  irw/ir0|  /,(£)- j^j/2(0)j 

(25) 

A*=™Mr0(( £)/,(/*)-(£)  /,(/})) 

(26) 

where 

(8-2 /Ao 

(27) 

and 

J|(£)~/  J}(x)e~**dx 

A) 

(28) 

WH-jgW) 

(29) 

(30) 

/,(£)  is  just  the  Laplace  transform  of  J2(x){4} 

/,</*)= £c,/2(»+^2) 

(31) 

where  Qx/2  is  the  Legendre  function  of  the 
kind  of  order  1  /2. 

is  therefore  given  by 

second 

(32) 

where  the  prime  indicates  differentiation  with  respect 
to  the  argument.  The  required  derivative  may  be 
found  from  the  recursion  relation{5) 

(X1  -  l)0i/2(*)  =  §  0,/2(*)-ie_,/2(*)  (33) 

For  convenience  in  evaluation,  both  £>,/2  and 
Q-i/2  may  be  expressed  in  terms  of  complete  elliptic 
integrals01: 

-[a^i)]'/’«[(^T)'/1|  (34) 


f  2  \'/2k\\ 

f  2  \ 

L+i)  *[ 

l  jt+i  / 

where  K(k)  and  E(k)  are,  respectively,  the  complete 
elliptic  integrals  of  the  first  and  second  kind  of 
modulus  k: 

K(k)—  f  /2(  1  —  A: 2 sin2  r )  1/2 dt  (36) 
A> 

£(*)=  r/2(l-*2sin2/),/2</f  (37) 

A) 

Values  of  K(k)  and  E(k)  may  be  obtained  from 
standard  tables  or  from  readily  available  computer 
software. 

L3(£)  may  likewise  be  reduced  to  an  expression 
involving  K(k)  and  E(k).  However,  for  most  practi¬ 
cal  cases,  the  factor  ( 8/r0 )2  by  which  /3(fi)  is  multi¬ 
plied  is  so  small  that  the  contribution  to  A/?  from  the 
term  proportional  to  f3(fi)  is  negligible. 

3.  RESULTS 

To  illustrate  the  changes  in  impedance  as  a 
function  of  the  liftoff  distance  /  and  the  conductivity 
o,  calculations  were  made  for  a  loop  of  radius  r0  =  1 .27 
cm  (diameter  of  1  in.)  at  distances  /  from  0.05  to  1.5 
cm,  and  for  condqctivities  from  0.1  to  4  times  that  of 
aluminum  (o0=3.8X107  mho/m).  These  results  are 
shown  in  Figs.  2  and  3  as  a  function  of  /  for  various 
constant  conductivities.  The  normalized  dimension¬ 
less  changes  in  impedance  AX/cjfir0  and  A/?/w^r0 
are  chosen  as  the  quantities  to  be  plotted.  For  all 
values  of  conductivity,  the  value  of  hX/upr0  is  seen 
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to  approach  a  large  negative  value  as  l  decreases, 
showing  the  known  decrease  in  total  inductance  as 
the  loop  approaches  the  plane.  As  /  becomes  large, 
&X/b>nr0  approaches  zero  as  required.  Similarly,  in 
Fig.  2  A/?/w/ir0  is  seen  to  give  a  large  positive 
contribution  for  small  /,  and  approaches  zero  as  / 
becomes  large. 

To  illustrate  the  effects  of  the  conductivity  on 
the  changes  in  impedance  for  several  constant  values 
of  liftoff,  the  results  for  the  same  loop  are  shown  in 
Figs.  4  and  5.  The  change  in  reactance  AX/upr^  is 
seen  to  be  very  nearly  independent  of  conductivity 
over  the  range  considered.  The  value  of  AA/wpr0, 
however,  is  seen  to  increase  for  lower  values  of  a. 
This  resistance  term,  of  course,  approaches  zero  as 
the  conductivity  approaches  that  of  a  perfect  conduc¬ 
tor. 

Both  the  variations  in  resistance  and  reactance 
can  be  combined  into  the  one  graph  shown  in  Fig.  6 
by  plotting  &X  versus  A  R.  The  solid  lines  thus  show 
the  change  in  impedance  as  the  liftoff  is  changed, 
while  the  dashed  lines  show  the  variation  with  chang¬ 
ing  conductivity  for  constant  liftoff  L  Figure  6  repre¬ 
sents  a  narrow  strip  of  a  conventional  normalized 
impedance  plane  showing  curves  of  constant  conduc- 


Fig.  3.  Change  in  normalized  reactance  versus  liftoff  distance. 


tivity  and  curves  of  constant  liftoff,  the  strip  corre¬ 
sponding  to  large  values  of  conductivity,  with  the 
scale  of  normalized  resistance  greatly  amplified  in 
relation  to  the  scale  of  normalized  reactance.  In  a 
complete  impedance  plane  representation,  the  curves 
of  constant  conductivity  and  curves  of  constant  liftoff 
would,  of  course,  converge  at  the  origin. 


0  12  3  4  cT/CJr 

Fig.  5.  Change  in  normalized  reactance  versus  conductivity. 
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Fig.  6.  Chang,  in  reactance  versus  change  in  resistance. 
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Table  I.  Comparison  of  Theoretical  Change  in  Inductance  to 
Mutual  Inductance  for  the  Perfect  Conductor  Case 


/ 

(cm) 

Ajf 

0*4/0 

AZ. 

<**H) 

M 

0*H) 

2.5 

1.0887 

0.01737 

0.01727 

5.0 

0.5425 

0.00*66 

0.00869 

15.0 

0.07784 

0.00124 

0.001 14 

The  limiting  values  of  &X/u>nrQ  for  large  values 
of  o  can  be  checked  by  comparing  the  calculated 
values  with  that  of  the  case  of  a  loop  above  a 
perfectly  conducting  plane.  Using  image  theory,  the 
mutual  inductance  between  two  identical  loops 
located  a  distance  2/  apart  can  be  calculated<6): 

A/=2.54iVr0 

where  N  is  a  tabulated  function  of  r0  and  /.  The 
values  of  M  and  AT  at  50  KHz  are  compared  in 
Table  I,  and  quite  good  agreement  is  found. 

4.  CONCLUSIONS 

For  the  commonly  occurring  case  where 
0.1  r0,  the  change  in  coil  inductance  is  essentially  the 
value  that  would  occur  if  the  substrate  were  perfectly 
conductive;  AT  is  thus  dominated  by  its  dependence 
on  liftoff.  The  change  in  resistance  is,  for  constant 
liftoff,  proportional  in  first  order  to  skin  depth  (or, 
for  constant  frequency,  proportional  to  the  square 
root  of  substrate  conductivity);  however.  A/?  is  also 
strongly  dependent  upon  liftoff.  Second -order 
changes  in  AT  and  A R,  due  to  small  variations  in  / 
and  o  about  nominal  values,  are  well  approximated 
by  linear  functions  of  A /  and  Ao;  hence  variations  in 
AT  and  A R  may  readily  be  interpreted  in  terms  of 
corresponding  variations  in  liftoff  and  conductivity 
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NONDESTRUCTIVE  EVALUATION  OF  DIELECTRIC  LAYERS  ON  CONDUCTIVE 

substrates  by  microwave  surface  electromagnetic  waves 

C.  0 .  Gardner,  S,  A.  Long,  and  W,  Ou 
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University  of  Houston 

Houston,  Texas  7700s 


Ab  s  l  rac_t 

The  thickness  .md  permittivitv  (or  dielectric  constant)  of  dielectric 
l.r.rrs  <n  elect ricallv  conductive  substrates  can  be  determined  by  suit¬ 
able  measurements  using  surface  el ect romagnet ic  waves.  The  approach  used 
here  is  t.i  measure  the  cutoff  frequence  of  the  TM]  SEW  mode  and,  in  effect, 
tiie  propag.it  ion  constant  of  the  TM]  SEW  mode  as  a  function  of  frequency. 
The  theory  of  the  method  and  some  preliminary  results  obtained  using 
8-IJ0H/  SEW  supported  bv  a  layer  of  po lypropy lene  on  an  aluminum  substrate 
aie  presented.  Bv  going  to  higher  frequencies  the  method  can  be  extended 
to  thin  protective  coat ings  on  metals,  e.g.,  ceramic  coatings  on  jet 
engine  and  coal  utilization  components. 


INI  K>M)t  e  1  ION 

The  ide  i  ij  using  stirfic.*  e  le«  t  romagnet  i  c  waves 
(SEWi  to  measure  the  thickness  and  dielectric, 
constant  for  tompley  permittivitv)  of  .»  I  aver 
of  dielectric  material  <*n  art  electrical.';  .in¬ 
ductive  substrate  is  well  known.  It  lias  re¬ 
ceived  i  considerable  attention  in  relation  to 
opt o - e l .  •  i  t ron  i  c  devices  in  wh i « h  thin  1  avers  of 
optically  transparent  material  on  conductive  or 
semi  conduct  ive  substrates  ate*  used  as  optical 
wa  vegu  i  des  .  ^  I  he  possibility  of  adapting 

the  technique  ter  nondest  i  uct  i  ve  evaluation  ot 
other  types  ot  <  oat  ings  wa  •  raised  bv  He'll  and 
coworkers.  However,  tin*  idea  does  not  appear 

to  have  been  pursued  to  t  fie  point  of  practical 
app I i 4  at  ion . 

There  are  currently  several  potential  applica¬ 
tions  of  the  method.  Several  in  particular  are 
noteworthy.  One  is  t  fie  case  of  protective 
coatings  on  components  of  jot  engines,  coal 
combustion  chambers,  magnet  oh vd rod vnamic  gen- 

(O 

orators,  and  the  like.  Another  is  the  case 

of  polymeric  coatings  for  environmental  protec¬ 


tion.  There  is  also  the  case  of  surfaces  of 
metals  prepared  for  adhesive  bonding,  where  the 
strength  attained  bv  the  bond  is  known  to  be 
sensitive  to  tilt*  condition  of  the  adherent  sur- 
f  a.  os . 

I iiis  paper  describes  some  preliminary  work  ex¬ 
ploring  some  of  the  possibilities  and  practical 
problems  associated  with  the  SEW  method.  It  is 
not  aimed  at  a  specific  application.  The  work 
involves  the  use  of  mic  rowaves  in  the  8-12  GHz 
range  (free  space  wavelengths  around  3  cm). 

"tii >at  ings"  are  simulated  by  relatively  thick 
(  1 .  on)  slabs  of  a  readily  available  plastic 
(polypropylene)  placed  on  a  large  sheet  of  alum¬ 
inum.  For  application  to  thin  coatings,  the 
techniques  would  have  to  be  "scaled  down"  by 
one  or  two  orders  of  magnitude.  Fortunately 
there  are  no  fundamental  obstacles  to  this. 

Thus  far  only  transverse  magnetic  (TM)  waves 
have  been  used;  the  possibility  of  using  trans¬ 
verse  electric  (TE)  waves  and  combined  TM  and 
TE  waves  remains  open. 
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Figure  1.  Dielectric  Layer  on  Metal  Substrate 
THEORY 


Referring  to  Figure  1,  the  illustrated  unbound¬ 
ed  planar  structure  can  be  shown  from  Maxwell’s 
equations  to  support  TM  plane  waves  propagating 
in  the  z-direction  (with  surfaces  of  constant 
phase  perpendicular  to  the  z-direction)  for 
which  the  z-component  of  the  amplitude  of  the 
electric  field  is  given  by  the  relations: 


r  EgSin  xx 

CKx_<£ 

Vx)  " 

1  *‘0<x-O 

**  Eq  (sin  x£)e 

£<  x 

(1) 

where 

Ko 

-  Y-k2)"2 

(2) 

K 


(<a2)1/2 


(3) 


6  is  the  propagation  constant,  and  k*  -  w^u0eO* 
c  is  the  permittivity  of  the  dielectric  layer 
of  thickness  1,  Conditions  for  such  a  surface 
wave  exist  when  the  relation 

k0  1  6  1  «l,\  (4) 
is  satisfied.  For  this  case,  both  x  and  x^  are 
real,  positive  quantities. 


The  remaining  components  of  the  elactric  field 
intensity,  as  well  as  the  y-component  of  the 
magnetic  field  intensity  follow  from  Maxwell's 
equations  (the  x- component  and  z-component  van¬ 
ish).  When  the  requisite  boundary  conditions 
on  the  fields  are  imposed,  the  following  dis¬ 


persion  relation  is  found: 

(xt)tan(xt)  -  (5) 

Equation  (5)  implicitly  defines  the  propagation 
constant,  6,  as  a  function  of  kQ,  £  and  c.  This 
transcendental  equation  cannot  be  solved  exactly 
in  closed  form.  If  Equation  (5)  is  rewritten 
in  the  form 

xtanx  -  £[(c-i)(k0«.)2-x2]1/2  (6) 

where 

X  «  X£ 

it  becomes  clear  that  the  solutions  of  the  equa¬ 
tion  correspond  to  the  intersection  of  the 
(circular)  curve 

y  -  el(e-l)(kQl)2-xY/2 

with  the  curve 

y  “  xtanx  . 

Multiple  solutions  to  the  equation  occur,  a  new 
branch  occurring  as  x  increases  by  a  multiple 
of  tt.  For  x  -  nu,  the  right  hand  side  of  Equa¬ 
tion  (6)  must  vanish,  ie., 


k 

n 


n* 


i(t-l) 


1/2 


n  -  0,1,... 


(7) 


The  values  of  k^  given  by  Equation  (7)  corre¬ 
spond  respectively  to  cutoff  frequencies 


Hi  m 

n 


n*c 


«U-1) 


1/2 


(8) 


-1/2 

where  c  -  (u^c^)  i*  the  speed  of  light  In 

vacuo.  The  fields  corresponding  to  the  solution 
of  Equation  (6)  for  n*<x< (n+1)*  are  described 
as  the  TM^  SEW  inodes.  The  cutoff  frequency  for 
the  TMq  mode  is,  of  course,  xero;  the  TM^  mode 
can  be  excited  for  any  frequency.  The  higher 
order  modes  can  be  excited  only  for  frequencies 
exceeding  their  respective  cutoff  f requencles . 

We  note  that  a  determination  of  the  cutoff  fre¬ 
quency  of  any  mode  except  TMn  determines  the 

1/2  U 
quantity  £(c-l)A 


If  the  left  and  right  hand  sides  of  Equation  (6) 
are  expanded  as  a  power  series  in  x  about  the 


value  x  •  n»t  an  approximate  solution  is  ob¬ 
tained  which  may  be  written  in  the  form 


(r/  *  i>(<^>*>2<vk„>2 


Equation  (9)  is  valid  for 

% 

k  >  k  . 

0  n 

Equation  (9)  may  be  written 


Ur-)2-U1/2  =■  K^UKk  -k  ) 

kq  tun 


Thus  if  the  propagation  constant  3  is  determined 

as  a  function  of  the  wave  number  k^  (for  fixed 

values  of  £  and  i) t  for  values  of  greater 

than,  but  near  k  ,  the  quantity  (t-l)f/c  may  be 
n  2 

determined  as  the  slope  of  a  graph  of  [  (fi/k  )*" 

1/2  U 

-1]  versus  k^;  the  graph  will  intersect  the 

k~-axis  at  the  value  k  ,  which,  in  turn,  deter- 

0  n  1/2 

mines  the  quantity  *(t-l)  .  Writing 


X+[1-4(%)Y'2 

sk  ,, 

2(— r)2 
n’t 


1  -  <&• 


Thus  £  and  j,  are  determined  as  functions  of  the 

experimentally  measurable  parameters  s  and  k  . 

n 

PRISM  METHOD  OF  LAUNCHING  AND  RECEIVING  SEW 

One  practical  means  of  launching  a  SEW  is  Otto's 
prism  method, ^  illustrated  in  Figure  2.  If 
the  angle  of  incidence  $  is  such  that  the  in¬ 
ternal  angle  of  incidence  0  exceeds  the  critical 
angle  for  the  priao-air  interface,  then  the  x- 
component  of  the  (now  complex)  propagation 
vector  for  the  field  below  the  prism  is  pure 
imaginary.  By  varying  the  angle  of  Incidence 
+  (and  consequently  the  internal  angle  of  Ind¬ 
uce  0),  the  ratio  B/k^  can  (for  appropriate 
of  kg)  be  made  to  assume  the  value 


Figure  2.  Prism  Coupling  Arrangement  for 
Launching  and  Receiving  Surface 
Electromagnetic  Waves 

necessary  for  a  SEW  on  the  substraLe.  At  i  ti  t  >, 
condition,  a  SEW  will  propagate.  Thus  ilio  con¬ 
dition  for  launching  a  PM  surface  wave  mu\  t«» 


f.1/2sint>  -  £- 
P  k. 


[i+(Y  o2u -k  )-y/2 

£  u  n 


where  £  is  the  relative  permittivity  of  the 
prism  material. 

If  the  wave  inside  the  prism,  and  incident  on 

the  prism-air  interface,  were  an  ideal  plane 

wave,  there  would  be  a  sharply  defined  internal 

angle  of  incidence  0,  at  which  a  TM  SEW  could 

n 

be  launched.  In  practice  the  incident  wave 
comprises  plane  waves  with  a  range  of  propaga¬ 
tion  directions  distributed  about  a  central  ray. 
Hence  as  0  is  varied,  the  amplitude  of  the 

launched  TM  wave  varies  and  is  maximum  fot  the 
n 

theoretical  value  of  0 .  By  measuring  the  value 
of  0  at  which  the  amplitude  of  the  TM  SEW  1m 
maximum,  as  a  function  of  k^  (or,  equi valent J y , 
the  frequency  of  the  incident  radiation),  B/k^ 
is  determined  as  a  function  of  k^,  and  Equation 
(10),  (12)  and  (13)  may  be  applied  to  determine 


EXPERIMENTAL  METHOD 

A  block  diagram  of  the  experimental  arrangem«*tu 
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Figure  3.  Block  Diagram  of  Experimental 
Arrangement 

is  shown  in  Figure  3.  A  sweep  oscillator  which 
has  a  frequency  range  from  8.0  GHz  to  12.4  GHz 
is  used  as  the  microwave  source.  For  launching 
a  surface  electromagnetic  wave,  the  prism  coupl¬ 
ing  technique  has  been  used.  The  45°  prisms 
were  made  of  paraffin  wax  having  a  measured 

permittivity  t  *  2.22;  the  base  size  of  the 
P 

prisms  is  20  cm  ;  20  cm.  In  order  to  minimize 
direct  pick-up  of  any  radiation  other  than  SEW 
modes,  two  microwave  absorbing  screens  were 
placed  behind  (or  ahead  of)  the  prisms.  A  gap 
of  9  cm  was  left  between  the  microwave  absorber 
and  the  specimen  on  which  SEW  propagated.  Each 
microwave  horn  has  an  aperture  8  cm  *  8  cm.  In 
order  to  get  a  far-field  pattern  we  have  to  set 

the  distance  between  the  transmission  horn  and 
2 

prism  r>2D  /A,  where  D  is  the  dimension  of  the 
aperture.  For  our  case,  r  has  to  be  greater 
than  45  cm.  A  divided  circular  quadrant  (not 
shown  in  the  figure)  with  a  radius  of  80  cm  was 
built  and  the  transmission  horn  mounted  on  it. 
For  the  receiving  horn  a  similar  scanning  device 
with  smaller  radius  (r  -  50  cm)  is  used.  The 
transmission  and  receiving  horns  can  be  Inde¬ 
pendently  scanned  through  90°, 

For  coupling  the  wave  in  the  prism  into  the  SEW, 
'here  must  be  an  air  gap  between  the  prism  and 
the  dielectric  layer.  Experiment  shows  the  op¬ 


timum  gap  for  the  best-coupling  is  h  *  A/2, 
where  A  is  the  free-space  wavelength  of  Lhe 
microwaves . 

An  aluminum  sheet  (alloy  f/6061)  of  size  8'  -8" 
is  used  as  the  conductive  substrate.  Several 
polypropylene  sheets  of  the  samt  size  but  dif¬ 
ferent  thicknesses  are  used  as  the  dielectric 
coating  material.  The  polypropylene  sheets  are 
layed  on  the  aluminum  sheet  and  clamped  in 
order  to  minimize  air  space  between  the  poly¬ 
propylene  and  the  aluminimum  sheet . 

There  are  two  basic  parameters  we  have  to  mea¬ 
sure,  namely  the  frequency  and  the  incident 
angle  in  the  air.  Before  making  quantitative 
measurements  we  have  to  scan  several  times  to 
determine  the  angular  range  within  which  the 
TMq  and  TM^  modes  propagate  with  maximum  ampli¬ 
tude  . 

In  the  actual  measurement,  as  the  external  an¬ 
gle  $  is  increased,  we  have  to  adjust  the  posi 
tion  of  the  prism  slightly  in  order  to  keep  the 
central  ray  of  the  incident  beam  near  the  edge 
of  the  prism  for  most  efficient  coupling.  Then 
the  coupling  angle  can  be  obtained  by  scanning 
the  incident  beam  from  +  ■  0*  to  4  m  80°  and 
measuring  the  angle  at  which  the  most  energy  In 
coupled  into  the  surface  mode.  The  microwave 
frequency  is  obtained  from  the  frequency  meter. 


Figure  4.  Surface  Wave  Intensity  Versus  Exter¬ 
nal  Angle  of  Incidence 
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RESULTS 

Figure  4  is  a  representative  graph  of  the  de¬ 
tected  signal  amplitude  as  a  function  of  the 
external  angle  of  incidence.  The  large  peak 
corresponds  to  the  TM^  mode;  the  smaller  second 
peak  corresponds  to  the  TM^  mode.  Table  1 
shows  an  example  ot  measurements  obtained  from 
a  specimen  ot  thickness  «  *  1.50  cm  and  the 
dielectric  constant  of  the  polypropylene  layer 
t  =  2.2 5.  There  are  two  modes  (TM^  and  TH^) 
which  can  be  propagated  on  this  structure.  The 
data  m  Table  1  is  the  result  of  the  TH^  mode, 
because  only  the  higher  mode  is  useful  to  de¬ 
termine  the  dielectric  constant  and  thickness 
of  the  laver.  The  table  basically  contains  the 
frequency  f  and  the  incident  angle  in  the  air, 
the  other  data  is  simply  determined  from  f 
and  $.  The  f roe-space  wave  number  is  k^  »  2nf/c, 
where  c  is  the  velocity  of  1  ight  in  free  space. 
The  incident  angle  in  the  prism,  t)  t  is  obtained 
from  the  incident  angle  in  the  air,  <J> .  There 
are  two  i.ises:  i)  For  i  <  n/4  then  0  =  n/4 

-sin  *  (sin  l*1 /4-$ ) /n  );  ii)  If  <p  >  n/4  then 
-1  r 

0  -  n/4+sin  (sin(|-n/4)/n  ),  where  n  is  the 

P  P- 

re fraction  index  ot  the  prisms  (n  *  /e  ). 

P  P 

The  most  important  parameter  we  must  know  is 
the  ratio  B/k^.  According  to  the  theory  of 
launching  SEW  by  the  prism  coupling  technique, 
the  incident  angle  in  the  prism  is  determined 
by  the  equation  b  =  knn  sinb;  therefore  if  we 

0  p 

know  0  and  n  ,  B/krt  can  be  obtained  from  the 
P  o 

equation  b/k_  *  n  sinb.  0/ k_  in  turn  determines 

0  P  2  1/2  0 

the  quantity  [  (3/kg)  -1] 

The  approximate  dispersion  relation  for  any  TM 

n 

mode  of  SEW  near  cutoff  is  given  by  Equation 

(9).  From  Equation  (10)  we  can  see  that  if  we 
2  1/2 

plot  (  (tf/kg)  -lj  versus  k^ ,  a  straight  line 
will  result,  and  the  slope  s  will  be  equal  to 
(e-l)/tt;  the  Intercept  will  be  equal  to  the 
cutoff  wave  number  k  .  A  representative  graph 
corresponding  to  the  data  in  Table  1  is.  shown 
in  Figure  5.  For  this  particular  case,  the 
lope  a  -  8.39  *  10_3m  and  kj  -  186.7  m*1. 
ir  •  and  k^  are  determined,  c  and  l  are  ob¬ 


tained  by  means  of  Equations  (12)  and  (13).  \ 

For  the  example,  we  have  €  -  2.16  and  (  -  1 .  ‘>6  j 

cm.  These  may  be  compared  with  the  value 
t  “  2.25  determined  by  the  standard  waveguid. 
method,  and  t  ■  1.50  cm  measured  with  a  pan  u! 
calipers . 


Table  1.  Representative  Data  for  TMj  SEW  on 

1.5  Inch  Polypropylene  Layer  on  Ali«- 
lnua  Substrate 


Figure  5.  Function  of  Surface  Wave  Propagation 
Constant  Versus  Wave  Nufcer  of  Ex¬ 
citing  Radiation 

CONCLUSION 

The  possibility  of  measuring  the  thlrknesu  and 
dielectric  constant  of  a  dielectric  layer  on  a 
conductive  substrate  by  measuring  the  propaga¬ 
tion  constant  of  a  TM  SEW  using  the  prism 
launching  method  has  been  demonstrated  In  a 
regime  appropriate  to  the  8-12  GHZ  frequency 
range.  To  handle  thinner  dielectric  layers  It 


2U 


will  be  necessary  Co  employ  Much  higher  fre¬ 
quencies  . 

A  number  of  important  points  remain  to  be  in¬ 
vestigated,  including:  (1)  the  effects  of 
pronounced  variations  in  the  thickness  of  the 
dielectric  layer;  (2)  the  effects  of  pronounced 
variations  in  the  dielectric  constant  of  the 
dielectric  layer;  and  (3)  the  effects  of  im¬ 
perfections  in  the  surface  of  th'*  conductive 
subs  t rate . 

In  order  to  rationally  optimize  the  experimen¬ 
tal  arrangement  it  will  be  necessary  to  develop 
a  detailed  mathematical  model  of  the  prism  SEW 
launching  arrangement. 

Finally,  it  would  be  worthwhile  to  investigate 
alternative  launching  arrangements  including 
gratings  (or  similar  periodic  structures)  and 
special  horns;  and  to  investigate  the  possi¬ 
bilities  of  using  transverse  electric  (TE) 
modes  as  well  as  both  TM  and  TE  modes. 
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THE  CHANGE  IN  IMPEDANCE  OF  A  SINGLE-TURN 
COIL  DUE  TO  A  VOID  IN  A  CONDUCTING  HALF  SPACE 
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Department  of  Electrical  Engineering 
University  of  Houstoif 
Houston,  Texas  77004 


Abstract 

The  problem  of  detection  and  characterization  of  a  flaw  in  a 
conducting  half-space  using  an  eddy-current  coil  oriented 
parallel  to  the  interface  is  examined.  An  expression  is  de¬ 
rived  for  a  first  order  approximation  for  the  change  in  com¬ 
plex  impedance  due  to  a  void  located  within  the  conducting 
medium.  The  overall  impedance  is  a  function  of  the  radius 
and  lift-off  distance  of  the  test  coil  and  the  conductivity 
of  the  material.  An  analytical  expression  is  derived  for  the 
change  in  impedance  as  a  function  of  the  electric  fields  at 
the  position  of  the  flaw.  It  is  found  to  be  an  integral  over 
the  volume  of  the  flaw  of  the  electric  fields  found  with  and 
without  the  flaw  being  present.  The  limiting  case  of  a  de¬ 
generate  point  flaw  may  be  examined  in  greater  detail  by 
allowing  the  field  in  the  presence  of  the  flaw  to  be  approxi¬ 
mated  by  the  unperturbed  field.  For  flaws  small  enough  that 
the  field  does  not  vary  much  over  its  volume  the  field  may  be 
even  further  approximated  by  using  just  the  value  of  the  field 
at  the  position  of  the  centroid  of  the  flaw.  Plots  are  shown 
to  illustrate  the  behavior  of  the  change  in  impedance  as  a 
function  of  the  radial  range  of  the  flaw  and  the  depth  of  the 
flaw  centroid,  using  previously  derived  expressions  for  the 
fields  for  the  unflawed  case. 


1.  INTRODUCTION 

In  previous  work,  an  approximate  ana¬ 
lytical  expression  for  the  change  in 
impedance  of  a  single  turn  loop  par¬ 
allel  to  and  near  a  conducting  half 
space  has  been  obtained  analytically 
for  certain  idealized  test  coil  and 
specimen  configurations.^  The  re¬ 
sults  of  this  previous  investigation 
have  been  used  to  calculate  the  change 
in  impedance  due  to  a  finite  flaw  in 
the  conducting  material.  An  analytical 
expression  has  been  derived  for  a  first 
order  approximation  for  the  change  in 
the  complex  impedance  due  to  a  small 


void  within  the  conducting  medium  as  a 
function  of  the  electric  field  at  the 
position  of  the  flaw.  A  knowledge  of 
the  fields  everywhere  in  the  conducting 
medium  due  to  the  impressed  source  cur¬ 
rent  of  the  loop  is  required.  The 
change  in  coil  impedance  AZ  due  to  the 
void  can  be  found  in  terms  of  an  inte¬ 
gral  of  the  field  quantities  £  and  2^ 
(the  fields  at  the  position  of  the  flaw 
with  and  without  the  flaw  being  pre¬ 
sent)  .  The  electric  field  inside  the 
conducting  material  can  be  obtained 
from  the  standard  boundary  value  prob¬ 
lem  approach  in  terms  of  an  integral 


over  a  separation  parameter.  A  series 
expansion  of  the  integrand  then  allows 
the  field  to  be  expressed  in  terms  of 
a  derivative  of  the  Legendre  function 
of  the  second  kind  of  half  order.  The 
field  expressions  substituted  in  the 
equation  for  AZ  give  the  complex  im¬ 
pedance  as  a  function  of  the  depth  and 
radial  position  of  the  flaw.  The 
radial  position  affects  only  in  the 
magnitude  but  changes  in  the  depth  re¬ 
sult  in  variations  in  both  the  magni¬ 
tude  and  phase.  This  can  be  explained 
by  the  fact  that  for  our  eddy  current 
approximations  the  radial  position  of 
the  loop  will  always  be  a  negligibly 
small  portion  of  the  wavelength  in  air. 
However,  the  depth  in  the  conducting 
material  is  not  small  compared  to  the 
skin  depth  and  therefore  affects  the 
phase  of  the  response. 

2.  BASIC  THEORY 

Consider  a  two  port  network  consisting 
of  a  transmitter  T;  a  receiver  R  and  a 
conductive  body  B  containing  an  inter¬ 
nal  void  of  volume  Vp  as  shown  in  Fig. 
1.  Then  identify  the  multiply  con¬ 
nected  region  bounded  by  the  surface 
SF  of  the  void;  the  closed  surface  ST 
surrounding  the  transmitter  and  the 
associated  guiding  structure  (a  portion 
of  which  coincides  with  a  standard  re¬ 
ference  plane  cutting  through  the  field 
guiding  structure) ;  a  similar  surface 
SR  for  the  receiver?  and  a  boundary  at 
infinity  S^.  Let  the  volume  enclosed 
by  this  surface  be  V,  and  let  £q  and 
Hq  denote  the  time  harmonic  fields  that 
would  exist  when  the  material  within 
V_  is  identical  to  that  within  Vn  and 
E  and  H  denote  the  field  actually  ex¬ 
isting  with  the  void  present. 

Since  the  volume  V  contains  no  sources 
and  the  material  within  V  is  identical 


with  or  without  the  presence  of  Vp,  the 
Lorent2  reciprocity  theorem  is  valid. 
With  the  bounded  volume  V,  £Q,  HQ  and 
£,  H  satisfy  the  same  set  of  equations, 


V*2  *  -jcoiigH  ?  V*H  *  ju>£  (r)£ 


By  the  Lorentz  reciprocity  theorem,  we 
have 


£*HQJ  *ds 


0 


where  ds  is  an  element  of  surface  which 
is  normal  to  and  directed  into  the 
bounded  volume.  The  surface  S  is  the 
sum  of  surfaces  which  include  S^,  S^, 

SR  and  SF 


But 


L  ■  L  *  L  ♦  L  *  I, 


Since  the  tangential  component  of  elec- 


M 


trie  field  intensity  vanishes  on  the 
metal  waveguide  by  the  boundary  condi- 

tion'  j boundary  surface  "  °»  and  8ince 
the  fields  satisfy  the  radiation  condi¬ 
tion,  the  integral  over  Sw  vanishes. 
Hence,  we  have 


IRoIT(Z12-Zi20> 


-J  (i0xfi-ixH0) -ds 


.  •  .  42, 


'ST+SR 


£xh0) *ds 


Given  the  assumption  that  the  trans¬ 
mitter  and  receiver  waveguides  operate 
in  their  respective  dominant  modes, 
currents  and  voltages  are  defined  such 


E  *  Ve;  H  *  Ih;  *  VQe;  HQ  =  XQ?i 


where  e  and  h  are  the  normalized  dom¬ 
inant  waveguide  modal  electric  and 


magnetic  fields;  that  is  j 
•ds  *  1,  we  have 


cross 

section 


(£0x5-Sx80).di  -  WW 


I  (iLxH-ExS  ) -ds  -  ID  v  -I_V_ 
js  u  u  Ro  R  R  Ro 

R 

If  Eq,  Hq  are  the  fields  with  a  current 
IR  impressed  at  the  reference  plane  and 
the  transmitter  open  circuited  (IT  *0) 
and  £,  H  are  the  fields  with  an  im¬ 
pressed  current  IT  and  the  receiver 
open-circuited  (IR-0) #  then 

-ITVT  +1.  V-  -f  (20x8-|x3  ).ds  -  0 

0  0  SF 

Now  «  Z12  IR  ,  where  Z^2  is  the 
transfer  impeSanfie  between  transmitter 
and  receiver  with  no  void;  and 

VR  -  zx2IT'  where  zi2  is  the  tran*fer 
impedance  with  the  void. 


^  ^12~^12  ^  *  *  I  I  { 

X  0  XR  XT 


(E0*fi-i!xH0)  *ds 


|  IEqx3-Jxho ] *ds  «  -J  li^xH  using  Gauss*  theorem 


-  T-V  f  v*  (£**3-£*Hn)dv 

ZR  T  'V  U  U 


R0  T  VF 


-  f- V’  I  IH-  (Vxg  )-30*  (V>U) 


R0  T  VF 

-  Hq • (7xi)+l- (VxHq) ]dv 


Rq  t  'vp 


IH*  (- jonJ0H0)  -Eq*(  ju)i.E) 


-  H0*(-juu0H)+E*  (jue0E0)  ]dv 


ju(E-£n; 


IroIt 


But  e  -  en  -  ^ 
0  o> 


^ -  f  £n*l  dv 

'v„  0 


E-en  =  -  12- 

0  U) 


Az. 


rVl  V*dv 

^R  T  >\1  V 


R0  T  VF 


where  £  and  are  respectively  the 
electric  fields  found  with  and  without 
the  flaw  being  present.  The  limiting 
case  of  a  degenerate  point  flaw  may  be 
examined  in  greater  detail  by  allowing 
the  field  in  the  presence  of  the  flaw 
(£)  to  be  approximated  by  the  unper¬ 
turbed  field  £q  over  the  volume  of  the 


*Z12“I7i^|v  lVr'z,j2dv 

R0  T  VF 

For  flaws  small  enough  that  the  fields 
do  not  vary  greatly  over  its  volume, 
the  expression  may  be  further  approxi¬ 
mated  by  using  just  the  value  of  the 
field  at  the  position  of  the  centroid 

of  the  flaw  (r  ,z  ). 

c  c 


1 


n 


AZ 


12 


°VF 

\X7 


lE0(rc'zc)r 


3.  FLAW  IN  A  HALF  SPACE 

The  basic  geometry  of  the  problem  is 
shown  in  Fig.  2  and  consists  of  a  single 
turn  coil  of  radius  "r^"  oriented  par- 
allel  to  and  at  a  distance  miH  above  a 
conducting  half  space  within  which  a 
flaw  of  volume  Vp  is  located  at  a  posi¬ 
tion  (rc,zc)  with  respect  to  the  origin 
of  the  co-ordinate  system. 


A .  (r , z)  —  i  is. 
♦  )q  2 


|  lB2(a)e-a*+C,(a)eazJ 


•  J^tarjda 


.  r 


a  Z 

C3(a)e  1  Jx(ar)da 


Using  the  four  boundary  condition  equa¬ 
tions  resulting  from  the  continuity  of 
2  and  H  at  z  -  1,  and  z  ■  0  the  four 
constants  can  be  found. 


Our  principal  interest  is  evaluating 
the  vector  potential  inside  the  con¬ 
ducting  half  space,  which  in  turn  gives 
the  fields  in  that  region.  The  simplest 

choice  is  to  evaluate  and  use  the 
III  J 

expression  for  A^  to  calculate  the 
time  varying  field. 


The  resulting  derivation  for  the  co¬ 
efficient  is 


Wl(ar0)e 


-at 


(1  * 


The  field  in  each  region  satisfies 
Maxwell's  equation  and  hence  the  wave 
equation  for  the  vector  potential  X. 

v2X  4  k2x  -  -u3iBpr.,sed 


Making  the  usual  quasistatic  approxi¬ 
mation  in  regions  I  and  II  (that  the 

k2X  term  is  negligible)  and  with 
2 

k  *  jujuo  in  region  III,  the  solution 
for  A  in  the  three  regions  is 


B1(a)e~azJ1(ar)da 


where  a  is  the  separation  variable  and 


q^  *  /a  +juiyo 


which  gives 


‘°ir°  L 


(ar) J^farg) 


e-alealz 

a+a, 


ada 


The  effective  range  of  the  integration 
variable  a  for  which  the  value  of  the 
integrand  is  significant  is  much 
smaller  than  the  quantity  wyo.  For 
such  cases  -  may  be  expanded  as  a 
power  series  *  in  » 


/u >uo 


q 

q+a. 


U+  H 

1  u.i 


(r,z) 


I 

I 

I 


1  7  ' 

/a  +ju>po 


✓i+  3Sp 


For  a  good  conductor  wyo>> a2  |— 1«  | 

°1 

-ft  .  ■■  «  f  I-  (ft_)  +  {ft— )  2-  /<*Y  ^  1 

a+al  aI  1  'a^  . J 

a^  ^  /jaiyo 

=  d+j)  J 

where  6  *  /2/ajyo 

tr-  V  <4°> 

<4)2  =  -  4  («c)2 


(£-)3  =  -  -*4^-  («o>)3 

“l  4 


The  vector  potential  in  the  conducting 
material  as  a  power  series  in  skin 
depth  is 


•J11"-'1'  •  >W6 


•  J1(ar0)e_a4e“lzadc. 

+  62|  i  Jx (ar) Jx (arQ) e~al 

aiz  2 , 

•  e  *  a  da 

"  fi3  \]  ^  Ji(or,Ji 

•  e'a*eal2  a3da  . 1 


(ar0) 


Using  the  Lorentz  condition,  $*A  • 

3  A 

-ue  -rz  t  for  the  present  case  with 
circular  symmetry,  v*A  «  0,  and  the 
scalar  potential  $  must  be  time  inde¬ 
pendent.  This  means  $  will  give  rise 
only  to  a  static  field  which  is  of 
little  interest.  Therefore, 

*  .  .  it  .  .)«», 


So  the  field  without  the  presence  of 
the  flaw  in  the  conducting  half  space  is 
given  by 

Unirr 

-  «  ; 

'  0 


*  U°^r°  a)  ea+j)z/4Ij  |  (l-j)6 


(ar)  Jx  (arQ)  e~afi*  ada 


e  a2da . J 


(ar) 


-  |  62J1(or) J1(arQ) 

«  V!o  ue(i+j)*/«(j  j~  6Ji 

,00 

•  (arQ) e“a£ada  +  j  63^  (ar) (arQ) 

,oo 

lada  -  J  62J1(ar)J1(ar())e 


-aJl 


-at 


•  a*da . ] 

This  series  of  terms  may  be  expressed  in 
the  form  of  a  Legendre  function  by  tak¬ 
ing  the  derivative  of  the  Laplace  trans¬ 
form  of  the  quantity  J1(ar)J1(arQ  )<2) 


’DO 

J1 

J0  A 


(ar) J. (arQ)e”alda  =  — - — 


t2+r2+r2 
°l/2  {  ^rrQ  ) 


f  (ar) J, 


(arQ) e~a^ada 

-all. 


d  l 


•  J1 (ar) Jl (arQ)e  da]  *  1(a) 


1(a) 


— — —  TF 
ir/Frj;  M  1/2 


2  2  2 
*  +r  +r* 

where  x  *  — - 


1(a) 


it(rr0) 


371  °l/2(x) 


From  the  recursion  relation  for  Legendre 
functions,  the  derivative  may  be  ex¬ 
pressed  in  terms  of  the  functions  them¬ 
selves. 


I 

I 


°l/2 (x) 


Thus 


XT7  [2  °l/2(x)_  2  °-l/2(x>1 


(x  -1) 


ICO 

J ^  (ar) (arQ) e”atada* 


Tr(rr0)3//2(xi-l) 


^2  ^1/2*X^“  1 


1(a) 


Neglecting  terms  of  order  6 *  and  above, 

lo»lV^U)e(i^)^ll+j]6l(0) 


For  a  small  flaw  located  at  (r  ,-z^)  as 

c  c 

shown  in  Fig.  2,  Eg  can  be  approximated 
as 


uoiro 


r(l+j)zC/<S[l+j]5I(a) 


Since  the  same  coil  is  used  as  the  pri¬ 
mary  coil  and  also  as  the  secondary 
coil  IT  »  IR  and 

ov 

AZ  -  if 


AZ 


VFuu0r^[I(a)  ]2e_22c/6e_32(zc/5'1t/4) 


From  the  expression  for  AZ  above  it  is 
evident  that  all  the  dependence  on 
radial  position  is  contained  in  the 
I {a)  term  which  contributes  only  to  the 
magnitude  of  the  change  in  impedance. 
The  depth  of  the  flaw  is  seen  to  affect 
both  the  magnitude  and  phase  of  the 
change  in  impedance. 

The  constant  phase  factor  comes 

from  the  first  two  terms  of  the  series 
expansion  of  X.,  This  factor  would 
shift  slightly  if  more  terms  were  in¬ 
cluded  in  the  series  expansion. 

4.  RESULTS 


is  plotted  versus  rc  in  Fig.  3  for  sev¬ 
eral  different  liftoff  distances  varying 
from  0.3  to  1  times  the  radius  of  the 
loop.  Each  of  these  curves  tends  to 
zero  as  the  flaw  approaches  the  point 
directly  under  the  center  of  the  loop  or 
as  the  distance  from  the  periphery  of 
the  loop  becomes  large.  Each  curve  at¬ 
tains  its  maximum  value  when  the  flaw  is 
directly  under  the  wire  of  the  loop  and 
the  value  of  this  maxima  is  a  function 
of  the  liftoff  distance.  This  behavior 
is  to  be  expected  since  the  induced  eddy 
current  distribution  under  the  loop  is 
zero  at  the  center  and  a  maximum  at  the 
position  of  the  current  carrying  coil. 
The  current  density  then  decays  as  the 
position  of  the  flaw  moves  further  away 
from  the  loop.  At  any  constant  radial 
position,  the  magnitude  of  the  eddy  cur¬ 
rent  density  becomes  larger  as  the  lift¬ 
off  becomes  smaller.  This  functional 
dependence  is  shown  in  Fig.  4  for  the 
case  of  rc  -  rQ. 

Since  the  axial  position  of  the  flaw 
affects  both  the  magnitude  and  phase  of 
the  change  in  impedance, the  real  and 
imaginary  parts  of  this  factor  have  been 
plotted  versus  flaw  depth  in  Fig.  5.  As 
the  flaw  becomes  more  deeply  imbedded 
inside  the  material  the  phase  can  change 
enough  so  that  the  reactive  part  of  the 
impedance  becomes  negative  near  zc/6  * 
it/4.  Similarly  the  resistive  part  passes 
through  zero  at  zc/6  *  tt/2,  but  the 
magnitude  has  decayed  to  an  almost  neg¬ 
ligible  value  for  flaws  at  this  depth. 

The  real  part  of  the  impedance  is  zero 
for  a  flaw  at  the  surface  (*c“0)  owing 
to  the  assumption  that  Vj,  is  infinitesi¬ 
mally  snail.  Tnus  when  Vp  is  located 
at  the  surface,  the  current  distribution 
is  not  affected. 


To  illustrate  the  change  in  impedance 
as  a  function  of  the  radial  location 
of  the  flaw,  the  function  [I(a))2r? 


For  constant  values  of  lift-off  and  depth 


of  the  flaw  the  magnitude  of  the  change 
in  impedance  is  shown  in  Fig,  6  for 
linear  scans  of  the  coil  near  the  flaw. 
Curves  are  shown  for  the  coil  moving 
parallel  to  the  y-axis  at  x*0,  x-O.Sr^, 
x~rQ,  and  x«1.5rQ  (see  Fig.  2).  The 
behavior  is  as  expected  with  peaks  at 
the  points  where  the  flaw  is  nearest 
the  current  carrying  wire. 
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